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A NEW MODEL TO ESTIMATE DAILY ENERGY EXPENDITURE 
FOR WINTERING WATERFOWL 

RICHARD A. MCKINNEY 1 ’ 3 AND SCOTT R. MCWILLIAMS 2 


ABSTRACT.—Current models to estimate daily energy expenditure (DEE) for free-living birds are limited to 
either those that use fixed thermoregulatory costs or those that more accurately estimate thermoregulatory costs, 
but require extensive and often logistically difficult measurements. Here, we propose a model based on basal 
metabolic rate (BMR), activity budgets, and site-specific energetic costs of thermoregulation that requires only 
simple measures of ambient temperature and wind speed to provide estimates of DEE. We use the model to 
calculate the DEE of Buffleheads ( Bueephala albeola) wintering at six habitats that afford differing degrees of 
protection from exposure within Narragansett Bay, Rhode Island. Bufflehead activity budget data collected during 
the winters of 2001-2002 and 2002-2003, along with average temperatures and wind speeds at the sites, were 
used to calculate DEE that ranged from 46.9 to 52.4 kJ/hr and increased with increasing wind speed. The 
energetic cost of thermoregulation composed as much as 28% of total DEE and increased with wind speed. Our 
DEE values were 13.4% higher, and thermoregulatory costs were up to 2X higher than those calculated using 
an existing model that incorporates fixed thermoregulatory costs. We also saw an increase in feeding activity 
with increasing wind speed; sensitivity analysis of the effects of wind speed and feeding activity showed that a 
1 m/sec increase in wind speed at our sites increased DEE by 2.5%, whereas a corresponding increase in feeding 
activity increased DEE by 4.5%. This suggests that in temperate winter habitats, increased feeding activity may 
have a greater impact on Bufflehead DEE than wind exposure. Site-specific model estimates of DEE could also 
provide additional insight into the relative contribution of environmental conditions and changes in waterfowl 
behavior to DEE. Received 27 May 2004, accepted 12 January 2005. 


The daily energy expenditure (DEE) of a 
species is the sum of basal metabolic rate 
(BMR), thermoregulatory requirements, and 
the energetic cost of daily activities such as 
feeding, locomotion, and social behaviors. 
Quantitative assessments of the daily activities 
of wintering waterfowl have been used both 
to identify important habitats for these species 
and to assess their response to changes in hab¬ 
itat quality (Fredrickson and Drobney 1979, 
Brodsky and Weatherhead 1985a, Baldassarre 
et al. 1988, Paulus 1988). Waterfowl activity 
budgets may be influenced by habitat type 
(Turnbull and Baldassarre 1987, Rave and 
Baldassarre 1989) and site characteristics such 
as food abundance, protection from exposure, 
and level of disturbance (Nilsson 1970, Jorde 
et al. 1984, Paulus 1984, Quinlan and Baldas¬ 
sarre 1984, Brodsky and Weatherhead 1985b, 
Miller 1985). Changes in waterfowl activity 
may also be tied to changes in DEE that result 
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from the influence of habitat characteristics. 
For example, increased exposure to cold and 
wind may increase thermoregulatory energy 
costs, and therefore require increased feeding 
to offset higher energetic costs (Bennett and 
Bolen 1978, Hickey and Titman 1983). Mod¬ 
els that allow comparison between the ener¬ 
getic costs of thermoregulation and specific 
waterfowl behaviors could be used to deter¬ 
mine the relative magnitude of these costs, 
and may also provide insight into the effects 
of habitat quality on the DEE of resident wa¬ 
terfowl. 

Traditional measures of DEE for birds from 
time-activity budgets use multiples of BMR 
to estimate energetic costs of activities, but 
may differ in how the thermoregulatory com¬ 
ponent of DEE is estimated (Weathers et al. 
1984). Early estimates of DEE included either 
a fixed cost of thermoregulation or one based 
solely on ambient temperature (Kendeigh 
1949, Schartz and Zimmerman 1971, Koplin 
et al. 1980). Models subsequently evolved to 
include a means to more accurately estimate 
thermoregulatory costs, but only by the exten¬ 
sive measurement of many variables (e.g., 
whole-body thermal resistance, forced-con¬ 
vective resistance), some of which may be lo¬ 
gistically difficult to obtain for free-living 
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wildlife (Pearson 1954, Stiles 1971, Walsberg 
1977). Weathers et al. (1984) proposed the use 
of standard operative temperature, or indices 
that allow single-number representations of 
complex thermal environments, to overcome 
some of these difficulties. However, while 
providing a much more rigorous estimate of 
thermoregulatory costs, this approach is lim¬ 
ited by the need for the construction and cal¬ 
ibration of taxidermic mounts, and may be 
best suited for aviary or well-controlled field 
applications. To date, researchers estimating 
DEE for free-living birds using published ac¬ 
tivity-based models are limited to either those 
that use fixed thermoregulatory costs or those 
that more accurately estimate thermoregula¬ 
tory costs, but at the expense of extensive and 
often logistically difficult measurements of 
many variables. 

Previous studies estimating DEE for win¬ 
tering waterfowl have employed models that 
use factorial increases of BMR and that as¬ 
sume a fixed cost of thermoregulation (Wool- 
| ey and Owen 1978, Albright et al. 1983, Mor- 
* ton et al. 1989, Parker and Holm 1990). For 
wintering waterfowl in northern areas exposed 
to low temperatures and high winds, thermo¬ 
regulation may compose as much as 80% of 
daily energetic costs (Ettinger and King 1980, 
Walsberg 1983). These costs may vary be¬ 
tween wintering habitats because of differing 
degrees of protection from exposure to wind 
and cold (Porter and Gates 1969, Goldstein 
1983, Bakken 1992). If estimates of DEE are 
! to be useful in assessing habitat quality for 
wintering waterfowl, they need to include 
some measure of the energetic cost of ther¬ 
moregulation based on local environmental 
conditions. 

Here, we present an activity-based model 
that includes habitat-specific measures of ther¬ 
moregulatory costs to estimate DEE of water- 
fowl in different habitats. Our model requires 
only simple measures of ambient temperature 
and wind speed, along with waterfowl activity 
budgets and morphological measurements. 
Thermoregulatory costs arc calculated by us¬ 
ing heat loss via conduction and convection 
as a function of temperature and wind speed 
to estimate the metabolic heat production re¬ 
quired to maintain body temperature (Birke- 
bak 1966, Goldstein 1983). Because of the 
ability to estimate site-specific DEE based on 


local conditions, the model may be useful in 
evaluating habitats that provide differing de¬ 
grees of protection from high winds and ex¬ 
treme temperatures. Model estimates could 
also be used to provide insight into the rela¬ 
tive contribution of environmental conditions 
and differences in waterfowl behavior to 
changes in DEE. 

In this study, we used our model to estimate 
the DEE of Buffleheads (Bucephala albeola) 
at six wintering habitats in Narragansett Bay, 
Rhode Island, that afford differing degrees of 
protection from exposure to wind and cold 
temperatures. Our specific objectives were to 
(1) compare estimates of DEE obtained using 
our model with those obtained using a previ¬ 
ously published model that incorporates a 
fixed cost of thermoregulation, and (2) ex¬ 
amine changes in DEE across the sites and 
determine the relative contribution of wind 
speed and waterfowl feeding behavior to 
changes in DEE. 

METHODS 

DEE site-specific thermoregulation mod¬ 
el .—Our model incorporating site-specific 
thermoregulatory costs into DEE for winter¬ 
ing Buffleheads (hereafter, SST model) con¬ 
sists of (1) a thermoregulatory component 
(EE Xhermoreg )—an estimate of the metabolic heat 
production required to balance heat loss from 
the bird to the environment through conduc¬ 
tion and convection, and (2) an activity com¬ 
ponent (EE Activjty )—an estimate of additional 
energetic costs resulting from specific daily 
activities of wintering Bufflehead expressed as 
multiples of basal metabolic rate (BMR). We 
sum these components to arrive at an esti¬ 
mated DEE. In our model, metabolic heat pro¬ 
duction includes resting energy expenditure in 
a thermoneutral environment (i.e., BMR) and 
the additional energy expenditure required to 
maintain thermal equilibrium. The model uses 
average temperatures and wind speeds that co¬ 
incide with activity budget sampling at the 
sites; DEE is reported in kJ/hr. 

Basal metabolic rates were estimated from 
those of 16 North American duck species 
summari/ed in McNab (2003). A plot of BMR 
versus body mass for these species gave the 
relation; BMR = 4.05M 07, \ where BMR is 
basal metabolic rate in ml 0 : /hr, and M is 
body mass in g. Estimates of BMR were con- 
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verted to kJ/hr using a conversion factor of 
18.8 kJ/L 0 2 , derived from the average com¬ 
position of the Bufflehead’s winter diet 
(Schmidt-Nielsen 1997). Body mass was ap¬ 
proximated at 450 g for males and 325 g for 
females (Gauthier 1993). 

Before calculating metabolic heat produc¬ 
tion, we first determined when this component 
of a Bufflehead’s DEE is necessary by com¬ 
paring ambient temperature with their lower 
critical temperature, or the temperature below 
which metabolic heat production is required 
to maintain body temperature (Schmidt-Niel¬ 
sen 1997). Lower critical temperature (LCT) 
was estimated by the empirical relation: LCT 
= 47.2M 0 18 , where LCT is in ° C, and M is 
body mass in g (Kendeigh 1977). We com¬ 
pared effective ambient temperature (T ef or 
the ambient temperature corrected for the ef¬ 
fect of wind speed; Siple and Passel 1945) to 
LCT to determine whether metabolic heat pro¬ 
duction would be required to maintain the 
duck’s body temperature. If T ef was less than 
LCT, we assumed that metabolic heat produc¬ 
tion was required to maintain body tempera¬ 
ture; we then calculated this energy require¬ 
ment and included it in the final DEE. On the 
other hand, if T ef was greater than the lower 
critical temperature, we did not include met¬ 
abolic heat production. Effective temperature 
was calculated using the relationship derived 
by Siple and Passel (1945): 

Tef = T h - (T h - T a ) 

X (0.474 + 0.239 X Vu - 0.023 X u), 

where T ef is the effective temperature (° C) 
used for comparison with the lower critical 
temperature, T b is body temperature (° C), T a 
is ambient temperature (° C), and u is wind 
speed (m/sec). 

If T ef was less than LCT, we used an em¬ 
pirical model to estimate metabolic heat pro¬ 
duction as a function of temperature and wind 
speed (Goldstein 1983): 

H T+lt = a + b\fu, 

where u is wind speed (m/sec) and H T+U is 
metabolic heat production (watts). The coef¬ 
ficient b is determined empirically from data 
summarized by Goldstein (1983) on seven 
species of birds (body size 13.5-3,860 g) by 
the relation: b = 0.0092M 066 X AT 032 , where 


M is body weight in g and AT is the difference 
between lower critical temperature and ambi¬ 
ent temperature in ° C. The coefficient a is de¬ 
termined under conditions of free convection 
(u = 0.06 m/sec) by the relation: 

a = H r — feV 0.06, 

where H T is an adjusted metabolic rate in 
watts at ambient temperature (Goldstein 
1983). We estimated H T using a heat transfer 
model proposed by Birkebak (1966) that cal¬ 
culates conductive heat loss from different an¬ 
atomical regions of the bird to the environ¬ 
ment using geometrical representations (e.g., 
head represented as a sphere, body represent¬ 
ed as a cylinder) and heat transfer theory (Ap¬ 
pendix; Birkebak 1966). Morphological mea¬ 
sures of body dimensions (Fig. 1) can be ob¬ 
tained from the literature (e.g., Belrose 1980, 
Gauthier 1993) or from measurements of mu¬ 
seum specimens. Average values for live Buf- 
fleheads (n = 4, obtained from the Connecti¬ 
cut Waterfowl Trust, Farmington, Connecti¬ 
cut) and Bufflehead study skins (n = 16, ob¬ 
tained from the Harvard Museum of 
Comparative Zoology, Cambridge, Massachu¬ 
setts) are summarized in the Appendix. Also 
summarized in the appendix are the equations 
drawn from Birkebak (1966), which were 
used to calculate metabolic heat production. 
For these equations, a heat transfer coefficient 
(k) of 0.102 cal/cm/° C was used for the entire 
body surface (Calder and King 1974). The 
thermal conductance of Common Eider ( So - 
materia mollissima) in water (i.e., wet thermal 
conductance) has been shown to be 57% 
greater than it is in the air (Jenssen et al. 
1989); therefore, we used a heat transfer co¬ 
efficient of 0.160 cal/cm/°C to calculate heat 
loss from the ventral body surface to the wa¬ 
ter. Metabolic heat production was calculated 
as: BMR + Q head + Q neck + Q b reast + Q body + 

Qventral surface, where BMR is baSal metabolic 

rate and Q is the heat loss term for each body 
component. 

Estimates of additional energetic costs re¬ 
sulting from specific daily activities (EE Activity ) 
were calculated by multiplying the proportion 
of time spent in a particular activity by the 
energetic cost of that activity. We used pre¬ 
viously reported multiples of BMR, summa¬ 
rized in Table 1, to calculate the energetic 
costs of activities by multiplying the propor- 
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FIG. 1. Body dimension measurement points required for input into the SST model to estimate DEE (see 
Appendix). A = head length, B = head height, C = head width. D = body width, F = body length, G = body 
height. H = neck length, I = neck width, J = neck height. 


TABLE 1. 
and fixed-cost 

Energetic costs as a multiple of basal metabolic rate (BMR) of activities used in the site-specific 
thermoregulation DEE models. 

Activity 

Operational definition 

Multiple of 
BMR 

Reference 

Dive 

Diving for food 

5.1 

dc Leeu w b)% 

Surface 

Surface and pause between dives 

5.1 

de Lecuw 1996 

Look 

Peering through the water at the cove bottom 

l.cS 

Woolc\ and Owen 1978 

Courtship 

Social display toward individual of the opposite gender 

2.4 

Albright et al. 1983 

Agonistic 

Hostile interaction between two individuals 

1.8 

Wooley and Owen 1978 

Swim 

Locomotion 

3.5 

Butler 2000 

lly 

Locomotion 

12.5 

Wooley and Owen 1978 

Preen 

Maintenance of feathers 

2.1 

Albright et al. 1983 

Alert 

Not moving, but actively observing surroundings 

1.8 

Wooley and Owen 1978 

Rest 

Not moving with bill tucked in feathers 

1.4 

Wooley and Owen 1978 
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tion of time spent in that activity by the cor¬ 
responding multiple of BMR. The contribu¬ 
tion of physical activity to DEE (Table 2; 
EE Activily ) was then calculated by summing the 
energetic costs of all activities in which Buf- 
fleheads engaged. 

DEE fixed-cost thermoregulation model .— 
Estimates of DEE were calculated using a 
method that incorporates a fixed cost of ther¬ 
moregulation (fixed-cost model; Morton et al. 
1989). In this model, the thermoregulatory 
component (EE xhermoreg ) is fixed and estimated 
at 5.9 kJ/hr (Morton et al. 1989). Additional 
energetic costs resulting from specific daily 
activities (EE Activity ) were calculated as in the 
SST model by multiplying the proportion of 
time spent in a particular activity by the en¬ 
ergetic cost of that activity. These two com¬ 
ponents were then summed to arrive at fixed- 
cost model estimates of DEE. 

DEE-habitat correlations. —We identified 
six Bufflehead wintering habitats within well- 
defined coves or embayments of the Narra- 
gansett Bay estuary. Included were two me- 
sotrophic, rocky- and sandy-bottom embay¬ 
ments (Sheffield Cove: 41° 29' 41" N, 71° 22' 
89" W; and Mackeral Cove: 41° 29'28" N, 
71°20'86"W), two mesotrophic soft-bottom 
coves (Coggeshal Cove: 41° 39'32" N, 71° 
20' 52" W; and Brush Neck Cove: 41° 41'47" 
N, 71° 24' 48" W), and two eutrophic soft-bot¬ 
tom coves (Apponaug Cove: 41° 41'40" N, 
71° 28'58" W; and Watchemoket Cove 41° 
48' 00" N, 71° 22' 75" W). Cove areas ranged 
from 18.6 to 86.1 ha, with an average of 42.2 
ha. Each cove supported consistent numbers 
of Buffleheads throughout the winter (Novem¬ 
ber through April); the median flock size at 
the six sites (determined by bimonthly cen¬ 
suses during the winters of 2001-2002 and 
2002-2003) was 18, ranging from 13 to 41. 
In winter, Buffleheads spend the majority of 
their time on the water and tend to favor shal¬ 
low water habitats (<3 m) in protected coves 
(Stott and Olson 1973, Gauthier 1993). They 
feed by diving to the cove bottom where they 
consume benthic invertebrates including crus¬ 
taceans, gastropods, and bivalves (Yocum and 
Keller 1961, Wiemeyer 1967, Gauthier 1993). 

We used focal animal sampling to quantify 
activities of Buffleheads at each of the study 
sites during the winters of 2001-2002 and 
2002-2003 (Altmann 1974). We completed 


965 observations on individual birds, resulting 
in over 80 hr of activity budget data. Obser¬ 
vations were randomly distributed over sam¬ 
ple sites and time during the daytime through¬ 
out the winter period when ducks were present 
(November-April). We chose individual 
ducks at random (i.e., observations began with 
the /th duck from the left in each flock, where 
i was a randomly generated number) and ob¬ 
served through a 32-60X spotting scope or 
through 10 X 50 binoculars for 5 min; behav¬ 
iors were categorized as dive, surface, look 
(i.e., peering through the water at the cove 
bottom), courtship, agonistic, swim, fly, preen, 
alert, and rest (Table 1). Preening included 
wing flapping, stretching, and scratching. 
Gender for each individual was identified 
when possible, except in rare instances when 
we were unable to distinguish between fe¬ 
males and first-year males that had not yet de¬ 
veloped breeding plumage (Carney 1992). 
Therefore, we report results for “males” 
(showing breeding plumage) and “females” 
(includes first year males). Activity data were 
collected using an observational software pro¬ 
gram installed on a laptop computer (JWatch- 
er, Animal Behaviour Laboratory, Macquarie 
University, Australia; http://www.jwatcher. 
ucla.edu/). Prior to analysis, data were aggre¬ 
gated into the following categories: feeding 
(dive, surface, look), social (courtship, ago¬ 
nistic), locomotion (swim, fly), maintenance 
(preen, alert), and resting. Each sampling 
event at a site consisted of 20-30 five-min 
observations; final data were averaged by 
sampling event and by site. 

Sensitivity analysis. —We used linear re¬ 
gression analysis of SST model estimates of 
DEE versus wind speed and feeding behavior, 
respectively, to assess the relative contribution 
of each to DEE. First, we estimated DEE us¬ 
ing average values of feeding activity across 
all sites, and plotted DEE versus wind speed 
over the range of wind speeds recorded during 
the study (i.e., feeding activity held constant, 
wind speed varied; regression equation: DEE 
= [1.1 X wind speed] + 42.1). Second, we 
estimated DEE using average wind speed and 
temperature across the sites and plotted DEE 
versus the proportion of feeding activity (i.e., 
wind speed held constant, feeding activity var¬ 
ied; regression equation: DEE = [43.1 X pro¬ 
portion of time spent feeding] + 17.3). In each 


TABLE 2. Thermoregulatory costs (EEn^^), energetic costs resulting from specific daily activities (EE Activity ), and daily energy expenditure (DEE) of Buffle- 
heads at six winter habitats in Narragansett Bay, Rhode Island, 2001-2003. Values were calculated using thermoregulatory costs based on site-specific temperature 
and wind speed (site-specific thermoregulation model), and a previous method using fixed thermoregulatory costs (fixed-cost thermoregulation model; Morton et al. 
1989). Values are in kJ/hr ± SD. 
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case, we used average values of temperature 
and all other activities in the model. Regres¬ 
sion equations generated from each analysis 
were used to estimate the relative contribution 
of wind speed and feeding behavior to DEE. 
For wind speed, we calculated the average 
percent increase in DEE per 1 m/sec increase 
in wind speed. For proportion of time spent 
feeding, we calculated the average percent in¬ 
crease in DEE per 5% increase in the propor¬ 
tion of time spent feeding. 

Statistical analyses. —Differences in the 
proportion of time spent on different activities 
by males versus females were investigated us¬ 
ing two-tailed Student’s Mests on data aver¬ 
aged across all sample sites. Site-specific time 
budgets were calculated by averaging individ¬ 
ual observations by sampling event and then 
by averaging sampling events by site. Propor¬ 
tions were arcsine-square-root transformed 
prior to regression analysis (Fowler et al. 
1998:87—88). Wind speed and temperature 
were averaged by sampling event and by site. 
Regression analysis and analysis of variance 
were used to assess the influence of environ¬ 
mental conditions on DEE and feeding behav¬ 
ior. Statistical analyses were performed using 
SAS (SAS Institute, Inc. 2001). 

RESULTS 

Estimates of DEE for wintering Buffleheads 
generated using the SST model averaged 49.0 
± 8.4 kJ/hr, or 1,176 ± 202 kJ/day, and dif¬ 
fered by up to 12% among sites (Table 2). The 
mean thermoregulatory component of DEE 
( EE Thermore g ; Table 2) was 11.7 ± 1.1 kJ/hr or 
23.9% of total DEE; EE Xhermoreg increased with 
increasing wind speed (r 2 = 0.61, P = 0.067). 
DEE did not differ between males and fe¬ 
males; however, thermoregulatory costs were 
higher for females (mean = 12.5 ± 1.2 versus 

10.9 ± 1.0 kJ/hr for males; t 5 = —7.2, P < 
0 . 001 ). 

The mean DEE (all sites) calculated using 
the SST model was 13.4% higher than that 
calculated using the fixed-cost model (Table 
2). The thermoregulatory component of DEE, 

5.9 kJ/hr, composed 13.7% of total DEE cal¬ 
culated with the fixed-cost model. 

Daily energy expenditures of Buffleheads 
calculated with the SST model increased with 
increasing wind speed for males (r 2 = 0.67, P 
= 0.046), females (r 2 = 0.64, P = 0.055), and 


54.0 A 
_ 52.0 
^ 50.0 
w 48.0 
q 46.0 
44.0 

42.0 4- 

0.0 1.0 2.0 3.0 4.0 5.0 

Wind speed (m/sec) 



0.0 1.0 2.0 3.0 4.0 5.0 

Wind speed (m/sec) 

FIG. 2. Correlation of wind speed with (A) DEE 
and (B) time spent feeding for Buffleheads (males and 
females combined) wintering at six coastal habitats in 
Narragansett Bay, Rhode Island, 2001-2003. Wind 
speeds are means of all sampling sessions conducted 
at a site. Error bars are ± SE. 

males and females combined (r 2 = 0.76, P = 
0.023; Fig. 2A). The proportion of time spent 
feeding by Buffleheads also increased with in¬ 
creasing wind speed (r 2 = 0.67, P = 0.047 
Fig. 2B). Estimates of DEE that were gener¬ 
ated using the fixed-cost model showed no re¬ 
lationship between DEE and wind speed. 

Buffleheads spent 75.7 ± 4.3% of their time 
feeding during daylight hours, and females fed 
more often than males (77.1 ± 5.4% versus 
74.2 ± 6.9%; / 545 = -2.6, P = 0.004; Table 
3). Males, however, spent more of their time 
engaged in courtship activities (2.39% versus 
0.43%; f 545 = 7.4, P < 0.001). Males and fe¬ 
males (combined) averaged 16.8% of their 
time engaged in locomotion and maintenance, 
and 4.5% of the time resting (Table 3). Buf¬ 
fleheads at Mackerel Cove spent the greatest 
proportion of time feeding and the least in all 
other activities, whereas those at Coggeshal 
Cove spent the least time feeding and the most 
in all other activities, except resting. Overall, 
Buffleheads spent between 0.3 and 2.6% of 






McKinney and McWilliams • WATERFOWL DEE MODEL 


51 


T3 

<U 

D. 

w 

T3 

C 

T3 

C 

c3 

u . 
1) 

3 >, 
e 3 d 

« q 

D. C 


<D 

5 « 

S 

£ :S 

« d 
W £ 
to rt 


T3 

o 

60 


E- 


_ i 

a 3 

in o 

+ q 


c o 
73 .c 
(U 

£ , 

£ 02 

rj 

Q- JJ 
* O 

« g 

£ =3 

p |° 

u 

73 

d ^ 

UJ C 

S 5 

£J 


ii 
»& 


X) s 

E & 
< £ 


q 

ri 


o 

d 


q 

d 


q 

ri 


m 



(N 



q 



q 



— 



r- 



rj 



d 



VO 






d 



Ov 



vC 



d 



o 

rf 

r- 

r» 


00 

Ov 

m 


VC 

in 

o 

CC 

vC 

O 

VC 

m 

c 

m 


rf 

00 

m 


q 

-H 

q 


q 

q 

q 

q 

q 

q 

o 

q 

q 

q 

oc 

q 

q 

q 

— 

d 

d 

d 

ri 

d 

o 

d 

d 

d 

ri 

ri 

d 

d 

d 

ri 

ri 

ri 

d 

d 

d 

+1 

+1 

+1 

+ 1 

+ 1 

+ 1 

+1 

+1 

+i 

+ 1 

+1 

+ 1 

+ 1 

+ 1 

+ 1 

+1 

+ 1 

+ 1 

+ 1 

il 

4-1 

in 

-t 

o 


-t 

r- 

rf 

ri 

cr 

r^ 


n 

VC 


in 

n 

oc 

in 

00 

m 

ri 

VC 

q 

q 

o 

q 

q 

q 

q 

q 

00 

o 

q 

q 

00 

q 

o 

q 

(N 

q 

o 

q 

d 

(N 

CO 

in 

d 

d 

d 

in 

d 

ri 

ri 

ri 

d 

d 

d 

d 

d 

VC 

d 

d 

d 


T i- 

Ov 

rl- 

n 

00 


(N 

o 

r- 

VC 


m 

rf 


n 

r- 

O 


oc 

n 

in 

Ov 

VO 

q 

q 

q 

q 

q 

q 

q 

rj 

rj 

q 

q 

q 

q 

00 

q 

q 

q 

q 

vO 

o 

CO 

d 

d 

d 

d 

d 

d 

ri 

d 

d 

d 

d 

d 

ri 

o 

d 

d 

d 

ri 

+ 1 

+1 

+1 

+ 1 

+ 1 

+1 

+1 

+ 1 

+ 1 

+1 

+ 1 

+1 

+1 

+ 1 

+ 1 

+1 

+ 1 

4-1 

+1 

+1 

+1 

o 

vO 

OV 

in 


00 

OV 

Ov 

Ov 

—i 

cr 

(N 

_ 

cr 

<N 

00 

cr 


VC 

rj- 

o 


q 

00 

(N 

— 

q 

00 

n- 

sC 

q 

q 

q 

q 

00 

ri¬ 

q 

q 

q 

q 

n 

rj 


d 

d 

Ov 

d 

oo’ 

ov 

OV 

Ov 


d 

d 

Ov 

d 

ve 

00 

VC 

d 

oc 

VC 

d 


vO 



00 





oc 

in 

r- 


r- 

Ov 


ov 

Ov 


r- 

in 

q 

00 

q 

(N 

q 

q 

o 

q 

rj 

in 

q 

Ov 

q 

q 

q 

q 

q 

00 

q 

q 

q 

d 

d 

ri 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

ri 

ri 

+1 

+ 1 

+1 

+ 1 

+1 

+ 1 

+ 1 

+ 1 

+ 1 

+1 

+1 

+1 

+ 1 

+ 1 

+1 

+ 1 

+1 

+ 1 

+ 1 

+1 

+ 1 

q 

Ov 


00 

00 

oc 

o 

q 

q 

n- 

r- 

in 

o 

Ov 

rf 

q 

00 

— 

-t 

— 

>n 

(N 

q 



r- 

d 

d 

ri 

d 


q 

q 

VC 

q 


ri 

q 

q 

d 

q 

q 


Ov 



Ov 





d 

VC 

d 

VC 

d 

VC 


d 

Ov 


od 

Ov 

m 


co 

O 

00 

rf 


rf 

Tf 

n- 

Ov 

n 

00 

n 

O 

sC 

VC 


00 

00 

CO 


— 

•—i 

— 

in 

00 

r- 


Ov 

r- 

n 

m 

q 

r^ 

q 

q 

— 

Ov 

q 

rl- 

00 

ri 

d 

d 

<N 

d 

d 

d 

o 

d 

d 

o 

d 

d 

o 

d 

d 

d 

d 

d 

o 

d 

+ 1 

+1 

+1 

+ 1 

+ 1 

+ 1 

+1 

+1! 

+1 

+1 

+1 

+1 

+ 1 

+1 

+ 1 

il 

+1 

+1 

il 

+1 

+1 

-t 

m 

o 

Ov 

r~ 

co 

m 

00 


cr 

00 

o 

00 

n- 

VC 

CT 

00 

vC 

OV 

cr 

— 

00 

in 

(N 

00 

q 

q 

o 


— 

n- 

— 

cr 

q 


OV 

q 

q 

rj 

q 

rf 

q 


o 

’ 1 

d 


ri 

d 

d 

ri 

o 

d 

d 


o 

d 

ri 

° 


ri 

o 

" 




q 



q 

q 

q 

q 


q 


q 


o 






q 

O' 

(N 

ri 

q 

q 

ri 


sC 

ri 

q 


q 


q 

d 

rj 

q 

Ov 

q 

CO 

in 

o 

d 


d 

oc 


— 



00 

— 

d 

— 

d 


d 

oc 

vC 

d 

d 

+ 1 

+1 

+ 1 

+ 1 

+1 

+ 1 

+1 

+1 

+1 

+1 

+1 

+1 

+ 1 

+1 

+ 1 

+1 

4-1 

+1 

+ 1 

+ 1 

4-1 

vO 

— 

q 

q 

q 

q 

q 

q 

q 

q 

— 

q 

00 

00 

q 

q 

rj 

q 

rj 

q 

q 

vC 

3 

d 

ri 

d 

ri 

00* 

00 

d 

d 

d 

d 

d 

ri 

oc 

oc 

d 

ri 

d 

d 

d 

vO 

oc 

r- 

r- 

p- 

r- 

r- 

VC 

r- 

oc 

00 

00 

r- 

00 

p- 

sC 

r~ 

r- 

r~ 

r- 

r- 


T3 
<D 

<u .£ 

o 5 -£ w -2 . 

*3 s £ 3 s £ 3 


U- U 


„ o 

u* o 


<u .3 

« 1 -£ 
£ £ 


£ 

S £ o S £ u S £ u 


0J ._ 


JU ^ 
I IL 




u 


> 

U 

Cu 

> 

U 

d 

> 

u 

o 

> 

u 

> 

y 

Ou 


o 

q 


< 

22 

u 


m 


> 

o 

H 


v .r 

« £ -g 


"3 

O 

£t 

G. 

s. 

"3 


2 £ u 


o 


\PI>CV Apponaug Cove. BRl'CV Brush Neck Cove: COC.CV Coggeshal Cove; MAKCV Mackerel Cove; SHFCV - Sheffield Cove; WATCV - Watchemoket Cove. 







52 


THE WILSON BULLETIN • Vol. 117, No. 1, March 2005 


their time in social activities, and this time 
decreased as the ducks increased feeding (r 2 
= 0.71, P = 0.043). Similarly, the amount of 
time spent in maintenance activities (range = 
3.7-9.7%) decreased as time spent feeding in¬ 
creased (r 2 = 0.96, P = 0.001). 

Sensitivity analysis of model estimates of 
DEE versus wind speed at constant proportion 
of time spent feeding showed that a 1 m/sec 
increase in wind speed resulted in a 2.5% in¬ 
crease in DEE. Analysis of DEE versus the 
proportion of time spent feeding at constant 
wind speed showed that a 0.05 increase in 
proportion of time spent feeding resulted in a 
4.5% increase in DEE. 

DISCUSSION 

Our estimates of DEE using the SST model 
for Buffleheads at the Narragansett Bay win¬ 
tering sites (1,175 ± 202 kJ/day) are higher 
than those predicted from the fixed-cost model 
(1,036 ± 202 kJ/day), which uses a single en¬ 
ergetic cost of thermoregulation. Thermoreg¬ 
ulatory costs predicted by the SST model con¬ 
stitute up to 28% of the animal’s total DEE 
and are approximately twice as high on av¬ 
erage as that used in the fixed-cost model. 
DEE estimates for Buffleheads at our sites 
were also higher than a field metabolic rate 
predicted by an allometric relation of energy 
expenditure based on empirical studies (606 
kJ/day, non-passerines; Nagy et al. 1999). 
However, many of the studies from which this 
relation was derived were carried out in the 
breeding season in warm ambient tempera¬ 
tures, so our higher DEE estimates may be 
attributed in part to environmental conditions 
and the inclusion of thermoregulatory energy 
costs. 

Our model does not include the contribu¬ 
tion of heat gained from solar radiation or heat 
lost through evaporative water loss because 
these effects are likely relatively small (less 
than 10% of heat loss; Scholander et al. 1950, 
Strunk 1971, Wolf and Walsberg 2000), and 
were likely similar between our study sites. 
Nonetheless, these constraints limit the appli¬ 
cation of our model to comparative, single¬ 
species studies between habitats that are lo¬ 
cated in a similar geographic region. It is also 
important to note that the SST model is lim¬ 
ited by the availability of empirically derived 
energetic equivalents of specific waterfowl be¬ 


haviors, as is the fixed-cost model. We applied 
the model to Buffleheads, but were restricted 
to using literature-based energetic equivalents 
that were not specific to that species. There¬ 
fore, the DEE estimates presented here, while 
higher than those calculated from the fixed- 
cost model and from body mass alone, fall 
well within the probable error of 20-40% pro¬ 
posed by Weathers et al. (1984) for models 
that rely on generic energetic equivalents. 
However, while it would be difficult to argue 
that our model estimates are more accurate 
than those calculated from fixed-cost or body 
mass models, the utility of our model lies in 
the ability to determine the relative contribu¬ 
tion of wind speed, temperature, and specific 
waterfowl behaviors to DEE across sites with 
different environmental conditions and levels 
of activity. 

Wintering waterfowl may incur substantial 
thermoregulatory costs depending on ambient 
temperatures and the combined effect of wind 
and cold, and these may lead to increases in 
DEE. Changes in the relative amounts of ac¬ 
tivities exhibited by wintering Buffleheads 
may also alter DEE. In our study, estimates of 
DEE calculated with the SST model were cor¬ 
related with wind speed (Fig. 2A). However, 
feeding activity also increased with increasing 
wind speed (Fig. 2B), which could also con¬ 
tribute to an increase in DEE. Sensitivity anal¬ 
ysis of the effects of increases in wind speed 
and feeding activity on DEE showed that in¬ 
creases in feeding activity resulted in a rela¬ 
tive increase in waterfowl DEE nearly twice 
that of a corresponding increase in wind 
speed. Feeding activity may increase because 
of decreased prey abundance, or because of 
changes in the availability or energetic content 
of prey. Further studies at our sites have 
shown that feeding activity increased with de¬ 
creasing prey abundance, and also with de¬ 
creasing prey energy density resulting from 
changes in available prey species at a site and 
inter-specific differences in the energetic con¬ 
tent of prey (RAM and SRM unpubl. data). 
However, other factors, such as intra- and in¬ 
ter-specific competition and increased ener¬ 
getic demands, may also influence the amount 
of feeding activity. Although we are uncertain 
as to the cause of increases in time spent feed¬ 
ing at our sites, our results show that increased 
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feeding activity may have a greater impact 
than wind exposure on DEE of Buffleheads. 

Increased feeding activity may also affect 
the short- and long-term survival of Buffle¬ 
heads. For example, if wintering Buffleheads 
need to spend more time feeding, time for oth¬ 
er activities such as courtship and pair for¬ 
mation may be limited (Drent and Daan 1980, 
Meijer and Drent 1999). Although they exhib¬ 
it long-term pair bond formation and a high 
degree of flock synchrony, which results in a 
relatively small proportion of time spent in so¬ 
cial behaviors, courtship and maintenance ac¬ 
tivities are still important for their overall re¬ 
productive success (Robertson et al. 1998). 
Our results indicate that as Buffleheads at our 
sites spent more time feeding, they had less 
time available for maintenance and social be¬ 
haviors, which may have an impact on both 
their short- and long-term survival. This, cou¬ 
pled with the greater increases in energetic 
costs due to feeding activity predicted from 
model sensitivity analysis, suggests that DEE 
of wintering waterfowl in harsh climates 
would be lower in habitats with both high 
prey density and adequate protection from ex¬ 
posure. For example, sites such as Brush Neck 
Cove, which had the highest prey abundance 
(RAM and SRM unpubl. data) and also the 
lowest thermoregulatory costs for Buffleheads 
(Table 2), may be better candidate sites for 
protection as waterfowl wintering habitats 
compared with sites such as Mackerel Cove, 
which had low prey abundance and high ther¬ 
moregulatory costs. 

In summary, our SST model estimated DEE 
as the sum of basal metabolic rate and site- 
specific energetic costs of activity and ther¬ 
moregulation. The primary benefits of the 
SST model compared to other approaches in¬ 
clude its ability to (I) evaluate the effect of 
thermoregulatory costs on DEE of wintering 
waterfowl using simple measurements of wind 
speed and ambient temperature, (2) predict the 
extent to which the behavior of waterfowl dur¬ 
ing winter affects DEE, and (3) track changes 
in DEE over different time scales (i.e., hourly, 
daily, or seasonally) if the corresponding ac¬ 
tivity and environmental data are available. 
Also, because of its ability to estimate site- 
spccilic DEE based on local conditions, the 
model may be useful in evaluating the quality 
of waterfowl habitats that have different attri¬ 


butes such as prey abundance, or degree of 
protection from high winds and extreme tem¬ 
peratures. However, further studies will be 
needed to establish the independence of be¬ 
havioral responses to environmental condi¬ 
tions from the primary effect of the conditions 
themselves on the DEE of resident waterfowl 
before model estimates can be used in habitat 
assessment. 
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APPENDIX. Variables used to calculate heat transfer, an adjusted metabolic rate at ambient temperature, 
using a heat transfer model proposed by Birkebak (1966). Representative values are from repeated measurements 
on live and preserved Buffleheads from northeastern estuaries. Equations are taken from Birkebak (1966); k is 
the heat transfer coefficient, AT is the difference between body temperature (39° C) and ambient temperature. 

Variable 

Symbol 

Equation 

Representative 
value (cm ± SD) 

Head length 

A 

— 

5.9 ± 0.4 

Head height 

B 

— 

5.0 ± 0.8 

Head width 

C 

— 

3.2 ± 0.4 

Body width 

D 

— 

9.1 ± 0.9 

Body length 

F 

— 

18.2 ±1.1 

Body height 

G 

— 

6.3 ± 0.6 

Neck length 

H 

— 

2.0 ± 0.3 

Neck width 

1 

— 

2.9 ± 0.4 

Neck height 

J 

— 

2.9 ± 0.4 

Integument depth-body 

AX h<H]y 

— 

0.4 ± 0.1 

Integument depth-head 

AX head 

— 

0.7 ± 0.2 

Integument depth-neck 

AX ncck 

— 

0.7 ± 0.2 

Inner radius of body 

l*i hotly 

6wv = (D + G)/4 

3.9 ± 0.5 

Inner radius of head 

l*i head 

>*, w = + 0/4 

2.1 ± 0.4 

Inner radius of neck 

neck 

r„ UTk = (/ + J)/2 

2.9 ± 0.6 

Length of body 

L'body 

L hodv = F - (D + G)/2 

10.5 ± 0.5 

Length of neck 

^ "'neck 

L ne , k - //-(/+ J)I2 

0.9 ± 0.1 

Outer radius of body 

body 

FHhhIy OWv A X l, od y 

4.6 ± 0.8 

Outer radius of head 

C head 

0 )hr<ul * the ad A X/,,.,,,/ 

2.8 ± 0.7 

Outer radius of neck 

l*o neck 

r ()nctk ~ r in,;k AX,,,.,* 

6.8 ±1.1 

Area of ventral surface 

A vs 

A,, I'bntiv ^ 2r ohf)(h 

96.0 ± 21.6 

Heal loss from head 

Qhc;til 

= (2tt X r, #w X r jW X k X A/V 


Heat loss from neck 

Qnctk 

( F )hrtut 1 ihctitl ) 

Qnra - <2it x L„, a X k X A /*)/( 1 n( r ()nri J r tHn *) | 

_ 

Heal loss from breast 

Qhrcasi 

Q„r,a„ (2tt X /-, IWv X r lM X k X A/V 

— 

Heat loss from body 

Qbtidy 

X L Wv X k X A /*)/( 1 n( r olHMh / /- |/vhA ) ) 


Heat loss from ventral 

Qv, 

CL k x A tl% x ( /; */;,) 


surface 

Heat loss from tail 

Q,ad 

Qw// Qhrrau 








